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A SCHEME FOR DERIVING DAY-TIME 


BOUNDARY-LAYER WIND PROFILES 


By F. B. SMITH and D. J. CARSON 


Summary. This scheme for deriving day-time boundary-layer wind profiles uses readily 
available parameters, namely: time of day, month, cloud amount, surface geostrophic wind and 
surface wind. The scheme allows for the effects of surface roughness, stability and diurnal and 
seasonal variations in boundary-layer depth. Typical profiles for both rural and urban sites are 
derived. 

The scheme is illustrated fully by example in the section on ‘Application’, which may be read 
independently of the theoretical background which is outlined in the previous section. An 
extension of the basic scheme to allow for thermal-wind effects is described. 


Introduction. This paper describes a scheme for estimating the mean wind- 
speed profile in a day-time atmospheric boundary layer in terms of the following 
parameters: time of day, month, cloud amount, the surface geostrophic wind 
speed and the mean 10-metre wind speed measured by a properly exposed 
anemometer. At times it may be necessary to replace the surface geostrophic 
wind by the gradient wind if the air trajectories have significant curvature or to 
adjust its value to allow for thermal wind effects which, on occasions, can 
be very important throughout the boundary layer. 

The scheme takes account of boundary-layer stability and surface roughness 
in determining the nature of the wind profile in the lowest few tens of metres 
and also allows for diurnal and seasonal variations in the depth of the boundary 
layer which affect the way in which the wind speed approaches the ‘free- 
stream’ value (the geostrophic- or gradient-wind value) at the top of the 
boundary layer. 

Two categories of profile are illustrated according to the surface roughness of 
the terrain over which the profile is to be estimated. One may be classed as a 
rural profile at the site of the anemometer set within typical English country- 
side (i.e. above a well-maintained level grass surface surrounded at a reasonable 
distance by hedges, trees, villages etc.); the other is for use over an urban or 
forested area. 

The second section gives the theoretical background to the basic scheme with 
reference to source material. The third section gives the procedure for applying 
the scheme and indicates its limitations. A worked example is also given and 
this section can be read without reference to the other sections. A further 
section describes how the thermal wind can be incorporated into the scheme 
when it is important. 
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Theoretical background 


Parameters needed to specify the profiles. Straightforward formulations of the 
wind profile throughout the whole of the boundary layer do not exist. The 
composite profiles presented here are obtained from a combination of the wind 
profile based on standard similarity-theory arguments applying well below the 
top of the boundary layer, the height of the boundary layer and the magnitude 
of the ‘free-stream’ wind at the top. 

The parameters needed to specify the whole profile are: 

(a) 2%, the surface roughness. The effective z, generally varies from one part 
of the wind profile to another. For the near-surface profile the very 
local nature of the surface (e.g. short grass) at the site is important but 
away from the ground Zz, represents an average roughness over an 
upwind area which is significantly affecting the wind profile at the 
height concerned. Since at greater heights the effective z, will have to 
take into account hedges, copses, buildings etc. it will generally be 
greater than the z, used closer to the ground. 

Four values of z) have been used to determine the final composite 
profiles. At the rural site, assumed to be a typical meteorological 
observing station, Z) is assumed 3 cm close to the surface, z) = 10 cm in 
a layer containing the anemometer at 10 metres and z = 30 cm above 
about 100 metres, as a typical z, for varied English countryside. For the 
urban profile, which is only applicable at heights well above the height of 
the roughness elements themselves, z= 1m is taken from about 
20 metres merging once more into the profiles for z) = 30 cm at heights 
in excess of 100 metres where the urban environment will no longer 
significantly affect the profile unless the city extends for more than a 
few kilometres upwind. 

Some measure of the degree of instability in the surface stress layer. 

The total incoming short-wave radiation, R. 

The depth, A, of the atmospheric boundary layer. 

The value of the 10-metre wind speed, u19, corresponding to a surface 

roughness Z = 10cm or, alternatively, the surface geostrophic wind 

(Vg) or gradient wind (Vgr). 
A stability parameter in common use in the surface stress layer is the surface 
Monin-Obukhov length, L, defined by 

os —uPecpT 
gkH 

where u, is the surface friction velocity, p is the mean density of the air, cp is the 
specific heat of air at constant pressure, T is the mean absolute temperature of 
the air, g is the acceleration due to gravity, k is the von Karman constant and 
H is the surface sensible heat flux which is assumed closely related to the 
incoming solar radiation, R. From his model for boundary-layer profiles in 
unstable conditions, Smith’ produced curves of u,/Vg as a function of the 
surface Rossby number, V,/fz) (where f is the Coriolis parameter), and a 
stability parameter, H/V,. Also, knowing u;z/u,as a function of z, L and Zz 
for the surface layers, Smith provided curves showing the interrelationships 
between R, Vg, u,, L, uyo(Z) and Z. By use of these relationships it is found 
that the parameter Vg/u,9(Z)) can be represented, to a high degree of approxima- 
tion, as a function of z) and L only and can therefore be used as a stability 
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parameter in place of L, for any given Z. For z = 10cm, Vg¢/t. = 2.0 
corresponds, in the absence of thermal-wind effects, to neutral stability, 1.75 
to L » —100m and 1°5 toL y —10m. 

A further result from Smith’s curves is that the wind speed u,9(Z)) at 10 
metres corresponding to a surface roughness Z, is proportional to the wind speed 
Uy) at 10 metres corresponding to a surface roughness of 10 cm, under the same 
external conditions of R and Vg. The constant of proportionality is, to a good 
approximation, independent of the Monin—Obukhov length, Z. In particular: 


U9(Zo) nt __ J} 0°86 for z = 0-3 m, 
ie (Zo) = 0°73 for Z = 1°0m. 
The new stability parameter Vg/u,9(Z9) is therefore readily replaced in the 
analysis by Vg/ujo. 
The presentation of the profiles is greatly simplified by these observations 
and the format 


uz 
ra = f(z, Zo Vg/uo, h) 
10 


is adopted, where Z, is a function of height reflecting the upwind variations in 
surface roughness, as discussed above. 


Determination of the stability parameter Vg/uyy. In some situations both uw, and 
V, will be known from anemometer readings and synoptic charts. The stability 
parameter Vg/u,) will then be known directly. In areas of complex topography 
this may be the only acceptable way of proceeding. Otherwise one or other of 
Uyg and Vg must be available and can be used together with the incoming 
solar radiation, R, to deduce the stability parameter. However, R is not 
generally readily available but can be estimated from three easily assessed 
parameters, namely the time of day, the month and the cloud amount. Figures 
1 and 2, based on four years of radiation data from Cambridge, show the 
dependence of R on the three parameters. A value of R for o—1 oktas of cloud is 
estimated from Figure 1 and is modified for the actual cloud amount from 
Figure 2. The derivation of these figures, first used by Smith et alii,* is based on 
information on the Cambridge radiation data discussed and displayed in Smith 
and Jeffrey.’ 

Figure 3 represents Vg/u,) as a function of R and wu, and is derived as a 
direct consequence of Smith’s! curves and the inferred relationship between 
Vg/uy9 and L. Alternatively Vg/u,) can be represented as a function of R and 
Vg by the directly related Figure 4. 

Figures 1, 2 and 3 (or 4) provide then the stability parameter Vg/u,) which 
will be used in the classification of the profiles. 


The wind profile in the surface stress layer. The wind profile in the surface stress 
layer is determined from the Monin-Obukhov similarity theory formulation, 
Ou ou 
Be = ke PmllL) 


where, following Businger and Dyer and Hicks,5 9m(z/Z) for neutral and 
unstable conditions is given by 


om(2/L) = (1-16 z/L)+, | 2/L| <1. 
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After integration this gives 
ig ty x(&z) —x(&z,0) 
Uyo(Z%) — X(E10) —x(Ez,0) 
uz 
si 1(Zo)ti0 





x(Ez) = In [ = ] 4+ gtan-! , 


&, = (1—16 2/L)*. 
The low-level profiles were determined in this way with V,/u,. used as the 
stability parameter for classification purposes. 
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FIGURE I——INCOMING SOLAR RADIATION, Ry, IN MILLIWATTS PER SQUARE CENTI- 
METRE FOR O—I OKTA OF CLOUD 
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FIGURE 2—INCOMING SOLAR RADIATION MODIFIED BY CLOUD COVER 


As deduced from unmodified value Ry from Figure 1. 











Vio at met. station 


FIGURE 3—CONTOURS OF V¢/u,) IN TERMS OF R AND ty 


Values of u;9 are in metres per second. 
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FIGURE 4—CONTOURS OF t;9/R IN TERMS OF R AND Vg 


The wind profile above the surface stress layer. To determine the wind profile 
from the top of the surface stress layer to the top of the boundary layer we use 
the two parameters V,/u,, and A and assume that the wind speed will attain 


the ‘free-stream’ value Vg at z = h. 

The value V,/u) is obtained either directly or by the procedure described 
above. 

To obtain a reliable estimate of hf is not an easy task. There are, generally 
speaking, no estimates of this depth readily available, although radiosonde 
ascents might be used with caution to provide an estimate at certain times. It 
should be noted that formulae based on the steady-state similarity approach 
applied to the whole layer are only of use for estimating h if the boundary 
layer is neutral and has been so for some time. Such formulae are quite in- 
adequate for use with the generally non-steady, diabatic boundary layer. This 
is because the depth is determined by mechanical mixing and the total amount 
of heat which has entered the layer, not just the heat received at any instant. 

The determination of h, in the absence of any direct information, has been 
somewhat simplified and reduced to the nomogram of Figure 5. Here h is 
determined as a function of time of day, month, cloud amount and the wind 
speed at 10 metres at the meteorological station site (i.e. 10 cm is the effective 
surface roughness for the profile at 10 m). The curves are based on the develop- 
ing boundary-layer model of Carson® which has been used with mean boundary- 
layer temperature profiles obtained by Hardy (private communication) from 
a study of Balthum ascents at Cardington over a period of 10 years. This 
type of approach was first developed for a practical scheme, designed and 
discussed by Smith,’ for estimating the vertical dispersion of a plume from a 
source near ground-level. Here we have extended the basic model to allow for a 


contribution to h due to mechanically, as opposed to thermally, generated 
turbulent mixing. 
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FIGURE 5—-NOMOGRAM FOR THE DETERMINATION OF THE DEPTH OF THE BOUNDARY 
LAYER, h 

















Knowing the height 4 at which uz/u,9 attains the derived value Vg/ujo it is 
necessary to carry a somewhat subjective extrapolation of the low-level profile 
for 7) = 30 cm to meet this upper-boundary condition. 

Finally a matching process has to take place between the profiles for the 
three values of z). The profile for z) = 3 cm is given most weight in the first 
few metres, appropriate to the surrounding grass-covered area. At 10 metres 
the actual anemometer wind, u 9, gives a fixed point and account is taken of the 
slope of the profile there corresponding to z= 10cm. Above about 100 
metres the composed profile for zy = 30 cm described in the previous paragraph 
is given greatest weight, appropriate to the effect of distant, more severe, 
roughness. Naturally there is a further element of subjectivity in this matching 
process but with so many fixed constraints (the value of z, at the ground, uo, h 
etc.) the expected errors arising from this are really quite small. 
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Over an urban area the matching process blends the profile for 73 = I m 
below 100 metres into the derived profile for 29 = 30 cm above about 100 
metres. Note that the urban profile cannot have meaning below about twice the 
typical building height. 

These profiles are presented in Figure 6 for several values of Vg/u,) and h. 





Near ‘neutral Slightly unstable Unstable 


Vo = 7 1-75 15 
Yio 
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0-6 0-8 1:0 12 1:4 16 18 20 08 10 1:2 1-4 146 18 06 08 10 12 1:4 16 
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FIGURE 6—EXAMPLES OF WIND PROFILES 





Below 100 metres pecked lines denote urban profile and continuous lines the profile at meteoro- 
logical stations. 


The application of the scheme. This section describes how to use the 
scheme without reference to the theoretical background given above. To 
clarify the process an actual example will be followed step by step; the example 
assumes the following conditions: 

(a) Time: 14 cmt; month: July; place: a typical meteorological observing 

station in central England, set in the country. 
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(b) 


(c) 
Step 1. 


Step 2. 


Cloud amount has been fairly steady through the morning and midday 
period at 4 oktas. 
Average wind speed at 10 metres (height of anemometer) is 4 m/s. 


Estimate the incoming solar radiation R. 

From Figure 1 read off Ry appropriate to o-1 oktas cloud for the 
month and time of day. 

In the example Ry ~ 68 mW/cm? for July, T = 14 cnr. 

Modify R, according to the cloud amount, using Figure 2, to give R. 
Example: Cloud is 4 oktas, Ry is 68, which gives R as 48 mW/cm?. 


Estimate the ratio of the geostrophic and 10-metre winds, Vg/u,9. If 
Uy) is known, then Figure 3 gives Vg/u,) in terms of R (found from 
Step 1) and wp. 
If uy) is unknown but an estimated value of Vg is obtainable from a 
synoptic chart, then Figure 4 may be used to give uy) and Vg/u,. If 
both u,) and Vg are known directly then this step may be omitted or 
used simply as a check. Clearly in difficult topographical situations 
where our scheme may be suspect, direct estimates of Vg and u,, may be 
preferred. 
In the example v4) = 4 m/s, 
From Figure 3 the value of Vg/u,) appropriate to R = 48, uy) = 4 is: 
Ve 
— w 16. 
Uy9 
With significant curvature of the air trajectory, Vg should be replaced 
by the gradient wind following standard procedures. 


Estimate the depth of the boundary layer, i.e. where the ‘free-stream’ 
wind is first achieved. 

If a significant synoptically induced inversion exists, this provides an 
upper limit to the boundary-layer depth. Figure 5 gives the depth in 
the absence of any such inversion or when the inversion is higher than 
that given by the Figure. 

Start with the time of day in the top right-hand box; move upwards to 
meet the appropriate month curve; move left into the top left-hand 
box to meet the appropriate cloud-amount curve; move down to the 
bottom left-hand box to meet the 4) curve and finally move to the 
right and read off the estimated value of boundary-layer depth, A. 
The ‘path’ taken for the Example is indicated by the dashed line and 
yields an estimated h = 1080 metres. 


Estimate the wind profile. 
To enter Figure 6, the two parameters Vg/u,) and h, found above, have 
to be known. 


y 


In the example —* — 1-6 and h = 1080 m. 
“10 

The profile has to be interpolated between the curves shown. 

In our specific example the profile lies between those given for 

Ve/uy = 1°5 and 1°75, and following the ‘tails’ up to h = 1000m 

(approximately). 
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The values of uz/u,) which may thus be estimated are shown in Table I. 


TABLE I-ESTIMATED VALUES OF Uz/tj9 
Ug/tyo 

Zz Velo = 15 «21°75 «16* — Uus® 

metres m/s 

I 0-78 o72 O76 3:0 

10 I I I 40 

40 rir T2906 «6 «| 46 

100 122 -1°93Ss«d1°27 , 

200 128 1°45) 1°35 ° 
500 1°42 1°64 1°52 

1000 1°50 1°75 


* Appropriate to the example 


The corresponding wind profile over the city can be inferred using the 
dashed lines on Figure 6. These should only be taken down to within 
two building-top heights above ground. Below that wind speeds are 
very variable and are strongly influenced by individual buildings or 
arrays of buildings. 

In the example the winds at 1 m and 10m cannot meaningfully be 
defined. At 40 m the average wind is expected to be 4:1 m/s, which is 
less than the corresponding rural-site speed by 0-5 m/s. At 100 m and 
above the speed is hardly affected unless the town extends for more 
than a few kilometres upwind. 


* The thermal wind. In previous sections any variation of the geostrophic 
wind with height through the boundary layer has been neglected. However, in 
many situations this is not justifiable largely because horizontal temperature 
gradients result in changes in the horizontal pressure gradient with height. 

Four simplified situations may be considered (see Figure 7). Situations (a) 
and (b) result primarily in a change of geostrophic wind direction with height 


“es gt 
Cold Warm 


Normal westerly situation Normal easterly situation 


FIGURE 7—-THE FOUR SAMPLE SITUATIONS 
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and we will not concern ourselves with these further. Situation (c) results in an 
increasing Vg with height whereas in situation (d), Vg decreases with height. 
The magnitude of these changes is given by 


g >| 
ST Ly Ip 
where strictly 97/dy is taken normal to the isobars but within a constant 
pressure surface. Usually this may be approximated by 07/0y in the horizontal, 
normal to the isobars with an error which is usually well under 1 m/s in 
[V_g(z = 1000 m)—V,(z = 0)]. 

Since the difference in Vg between the surface and the top of the boundary 
layer is of primary importance, a representative 07T/0y appropriate to the whole 
depth of the boundary layer has to be estimated. Unfortunately screen tempera- 
tures (at approximately 1-5 m) are not very representative of the mean tempera- 
ture in the whole layer. Consequently some adjustment has to be made 
allowing for heat flux and wind speed and by using the conclusions of similarity 
theory which give typical temperature profiles as a function of stability and 
height. These factors are employed to give the potential temperature difference 
between the screen and z = 100 m in terms of R (the incoming solar radiation) 
and the wind at 10 metres. It will be assumed that the implied temperature at 
100 metres yields the required representative temperature of the layer. 

Probably an adequate procedure to estimate the thermal wind in the absence 
of radiosonde data is then to select two well-sited stations lying across wind and 
one or two hundred kilometres apart, take their screen temperatures, adjust 
them to 100 metres (using Figure 8) and take the station difference. 


‘ie / 














| | 
20 30 40 50 


R (mW/cm2) 





FIGURE 8— CONTOURS OF 0(1+5)—6(100) 
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If we denote the temperatures at the two stations by suffixes A and B, and let 
the stations be n hundred kilometres apart, then define 


_ Ta—Ta + 304— 205 , 





Axo00 T n 


where 
80 = 6(1°5) —0(100). 
The difference in the geostrophic wind between 1000 metres and the surface is 
then 
AVg = 3°3Aj0007 in United Kingdom latitudes. 
Given that the depth of the mixing layer is not 1000 metres but h metres the 
change in geostrophic wind speed through the layer is simply 


AnVg = = 3°34n7, 


1000 
defining A, 7. The magnitude of AnT is given by Figure 9. The conclusion of 
this simplified approach is that the magnitude of the geostrophic wind at the 
top of the mixing layer, /, in terms of uy9 is given in terms of stability (as in the 
earlier sections) modified by an additional increment AyVg, 

ie. Vg(z = h) = Vg(z = 1 m) + AnVy, 
remembering that V,g(z = 1 m) is a function of u49, Zz and R. The effect of 
AnV, is therefore to modify the ratio Vg/uj) used on page 249. Figure 10 
shows how this modified ratio can be found very easily. It is immediately 
obvious that Vg/u,) can now lie outside the comparatively narrow range, 1°5 to 2, 
determined by stability alone. 
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Figure 11 shows two sample profiles when the thermal wind effect is quite 
marked; the one on the left corresponds to an enhanced Vg due to cold air to the 
left of the airstream whilst the one on the right shows a nearly constant wind 
with height arising from warm air to the left of the airstream. 

Plotted in this way, profiles are readily constructed. The following constraints 
and steps may be used: 

(a) uz/uy = 1 at z= 10m. 

(b) At z= 1 mover grass u,/u,9- 0°75 with a gradient which if extrapolated 

up to z = 10 m would give a velocity ratio of 1-4. 

(c) The upper part of the profile is first approximated by a straight line 
that passes through 0-86 at 10 metres and Vg/u,) at z = 1000 m. This 
is then modified by smoothly running from the straight line into Vg/u,9 
at z = h, with 0u/dz = o there, as shown in the several examples. 














1 = ee oe ee ee A ee ee ze 
06 08 10 12 14 16 18 20 2-2 2-426 28 3006 08 10 1-2 
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FIGURE II-—-TWO EXAMPLES OF WIND PROFILES STRONGLY AFFECTED BY THERMAL 
WINDS 


Both examples are profiles at meteorological stations. 
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(d) Finally a smooth profile is constructed linking the parts and subject to 
all the constraints in (a) to (c). 

This rather bold approach must inevitably involve errors, but since the 
constraints are numerous, these must be relatively quite small. The effect of 
stability on these ‘mixed’ profiles therefore must be most pronounced in 
determining V,/u,), and the shape otherwise is not seriously affected. 
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A MEMORABLE RAINFALL EVENT OVER SOUTHERN ENGLAND 
(PART I) 


By PAULINE M. SALTER and C. J. RICHARDS 


Summary. During the period 14-16 September 1968 parts of southern England experienced 
a memorably heavy fall of rain and substantial flooding occurred in places; the rain system was 
particularly intense over the south-east during the 15th. From a meteorological analysis of the 
situation it appears that the dominant mechanism governing the precipitation over the south- 
east was a pronounced convergence zone in the lower troposphere, overlaid at upper levels by a 
strong positive vorticity advection field operating within a narrow zone of marked horizontal 
wind shear. This system became slow moving on the 15th and, within it, a low-level inflow of 
potentially unstable, moisture-laden air is envisaged as leading to the ‘dumping’ of a heavy fall 
of rain over a relatively small area. 


Introduction. Between 14 and 16 September 1968 parts of southern England 
experienced an unusually heavy and prolonged rainfall that resulted in extensive 
flooding in the south-east. Headlines appearing in the national Press on the 
16th, following the event, were vivid in their description and may be recalled. 
‘The Great Deluge—Garden of England is a Paddyfield’ (Daily Mail), ‘Floods 
Devastate South—Hundreds of homes evacuated’ (Daily Telegraph), “September 
in the Rain’ (Daily Express). In papers by Bleasdale! * this rainfall has been 
compared with other exceptional rainfall events in the United Kingdom. The 
purpose of this paper is to describe the sequence of events and the rainfall 
distribution and, from available data, to offer some explanation for the weather 
experienced in the light of the meteorological situation prevailing at the time. 
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Surface analysis and weather sequence. ‘Two rain systems affected 
southern England during the period 14-16 September and it was the second of 
these systems that precipitated the memorably heavy fall over south-east 
England; both systems were associated with a slow-moving depression in the 
vicinity of the Bay of Biscay. The first rain belt, characterized by vigorous 
thundery activity, moved northwards across south-west England during the 
morning of the 14th, reaching the London area by afternoon. The second, and 
major rain-producing system, developed over northern France and the English 
Channel during the evening of the 14th and moved slowly north towards the 
Channel coast of south-east England on the 15th. 

Weather over much of the British Isles had been dominated during the first 
half of the preceding week by a large area of high pressure over Scandinavia. By 
the 13th a rather slack pressure gradient was covering the U.K. A slow fall in 
pressure set in over and to the south-west of the British Isles during the 13th 
and accelerated over the south-west during the evening, and a depression with 
central pressure of a little under 1000 mb formed off Brest by midnight on the 
14th (see Figure 1(a)). Thundery rain and thunderstorms broke out over Brittany 
and the South West Approaches accompanying this development, and extended 
slowly north-eastwards over Cornwall and Devon in association with a cold 
occlusion. A noticeable feature of this system was the vigorous electrical 
activity associated with it, as registered in the form of atmospheric reports 
(SFLOCS). During the morning of the 14th the trough moved northwards 
over southern England, storms on and ahead of the trough triggering off very 
high repetition rates of SFLOCS in places, and by midday the thundery 
activity had transferred into the south Midlands, reaching south-east England 
and East Anglia by afternoon. A clearer, rather unstable south-westerly 
airflow behind the trough brought drier conditions into extreme southern 
counties during the morning. During the afternoon the trough became slow 
moving from the Bristol Channel to the Thames Estuary, and thundery activity 
started to die out. 

The development of the second system was heralded by a general fall in 
pressure over central France during the afternoon. This extended northwards 
in the evening, sharpening and rejuvenating the broad trough of low pressure 
then lying over the English Channel. The surface wind backed and strengthened 
north-easterly over the British Isles, promoting the development of a marked 
trough and zone of convergence over the Channel east of Cherbourg later on the 
14th (Figure 1(b)). Renewed rain and thunderstorms broke out over France 
during the afternoon, drifted slowly north over the Channel in the evening and 
reached Portland Bill and the Isle of Wight shortly after 00 cmr on the 15th. 
The rain subsequently spread over central southern England and south-east 
England, south of the Thames. By 06 cmt on the 15th the depression had 
drifted southwards over Biscay, extending a broad trough north-eastwards over 
the Cherbourg peninsula and the Dover Straits. The axis of this trough, 
marking the main convergence zone, transferred slowly northwards over 
central Kent, and remained slow moving in a north-east to south-west direction 
during the morning and afternoon, the strong north-easterly airstream to the 
north of the trough contrasting markedly with the slack southerly flow to the 
south (Figure 2(a)). Lying south-east of the trough axis, the extreme east of Kent 
and Sussex escaped the prolonged heavy rain that was affecting places to the 
north and west. ‘Though thunderstorms were quite extensive within the 








a 


i 


ae & oid 
SiC LbG thy trag 


2 


Pn — 
WT 








L 





FIGURE I—SURFACE ANALYSIS AND I000—-500-MILLIBAR THICKNESS, 14-16 
SEPTEMBER 1968 
(a) oo cmt, 14 September 1968. Thickness values are in geopotential decametres. 
Hatching shows wet air at 700 mb (dew-point depression < 2 degC). Stippling shows dry air 
at 700 mb (dew-point depression > 20 degC). 


convergence zone during the morning, the electrical activity (as observed in 
SFLOC reports) was by no means as intense as it had been over central 
southern England the previous day. 

A temporary improvement was occurring south of the Thames by early 
afternoon of the 15th, but renewed heavy rain broke out over north Kent, 
E.sex, and Suffolk during the afternoon as a small depression formed within the 
convergence zone over the Thames Estuary. This depression drifted south-west- 
wards over west Kent, returning heavy rain to parts of Kent and Surrey later in 
the afternoon (Figure 2(b)), but it disappeared as a separate entity near the south 
coast during the evening. The rain area extended slowly northwards over East 
Anglia, precipitating continuous moderate or heavy rain here till 03 Gmr on the 
16th. It also moved over the south Midlands but rain here was much weaker in 
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FIGURE 1—continued 


(b) 00 cmt, 15 September 1968. Thickness values are in geopotential decametres. 
See also notes under Figure 1 (a) 


intensity compared with that nearer the trough. During the evening of the 15th 
the axis of the trough drifted slowly westwards over the south-east allowing the 
clearer, drier conditions behind it to spread over Surrey, Kent and Sussex. A 
further small depression formed in the trough off the Norfolk coast by midnight 
on the 16th, and drifted over Norfolk, disappearing south of Cambridge by 
07 GMT, but at the same time, displacing the main zone of convergence west- 
wards over East Anglia, which it weakened considerably. By midday on the 
16th the depression over Biscay had moved over France and degenerated into 
a shallow, low complex system. 


Rainfall distribution. On 14 September rainfall totals in excess of 10 mm 
were recorded over the area of south-east and central southern England south of 
a line from Cheltenham to Colchester. The heaviest totals recorded on this 
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FIGURE 1—continued 


(c) 12 GMT, 15 September 1968. Thickness values are in geopotential decametres. 
See also notes under Figure 1(a). 


‘rainfall day’ (in the 24 hours from 09 Gt on the 14th to 0g GMT on 15th) were 
in a belt running north-eastwards from the east Hampshire border to just east of 
Maidstone in Kent. Several point falls greater than 100 mm were recorded 
near Tonbridge and Maidstone. A parallel belt of heavy falls from Croydon to 
Maldon in Essex had point falls over g0 mm but none in excess of 100 mm. 

Rainfall totals in excess of 10 mm covered a wider area on 15 September ;* 
most of the area south-east of a line from Weymouth to Spalding, including the 
whole of East Anglia, had more than 25mm. Again the belts of heaviest 
rainfall ran in a south-west to north-east direction. The widest belt with totals 
greater than 50mm, from Winchester to the Thames Estuary, including 
Clacton in Essex and Herne Bay in Kent, was approximately 50 km wide and 
had seven centres with recorded rainfall totals of more than 100 mm. These 





*In this section dates refer to the ‘ rainfall day’. 
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FIGURE 1—continued 


(d) 00 cmt, 16 September 1968. Thickness values are in geopotential decametres. 
See also notes under Figure 1(a). 


included the highest totals recorded for the rainfall day, 129.5 mm at Bromley 
in Kent and 123.7 mm and 124.5 mm at Godstone in Surrey. The second belt 
with totals greater than 50 mm was narrower, 15-25 km, over most of its length 
from Beaconsfield in Buckinghamshire widening to about 50 km as it crossed 
Norfolk and Suffolk to the coast with Great Yarmouth on the axis. Although a 
large part of this belt had recordings of more than 75 mm, there were only five 
reports of more than 100 mm. 

The combined distribution for these two days from 09 Gmt on 14 September 
to 0g GMT on 16 September (Figure 3) shows the area of heaviest falls from 
Hampshire, across north Sussex, Surrey and Kent and across the Thames to 
Essex. The highest two-day totals were reported in Essex—201.4 mm at Marsh 
Farm Sewage Works and 200.9 mm at Stifford Pumping Station—in the 
Tilbury—Purfleet area. (Other reports are given in Table I.) In these two days 
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FIGURE 2—continued 
Shading shows the rain area. 
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FIGURE 3—TOTAL RAINFALL FOR 14-15 SEPTEMBER 1968 
National Grid Lines are shown at 100-km intervals. 


most of south-east England and East Anglia had more than 5 per cent of its 
average annual rainfall. Part of Kent and Surrey, about 450 km?, and an area 
of about 750 km? along the Thames Estuary had more than 20 per cent of their 
average annual rainfall. The highest percentages for the two days were in 
Essex at the two stations recording the highest rainfall totals with 39-0 and 39°7 
per cent respectively. A further small amount of rain fell in the area on 16 
September but totals were mainly less than 10 mm and were almost negligible 
compared with the amounts on the previous two days. 

Autographic records show few intense short-period falls during the two days 
but do show that most of the two-day total of rain fell in 24 hours or less, 
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spanning the two rainfall days. The Broadfield, Crawley, chart shows the main 
rainfall beginning between 05 and 06 crt on the 15th and continuing steadily 
until 13 mT followed by slight rain for about 5 hours. At 18 GMT it began 
to rain steadily again and this lasted for 3 hours. This was a fairly typical 
pattern for stations in Kent and Surrey. North of the Thames the chart for 
Basildon shows two storms starting at 1630 GMT and 22 cmt on the 14th and this 
was followed by periods of slight rain until the end of the rainfall day. On the 
15th it rained steadily for 5} hours from 10 cmt with further outbreaks of rain 
after 20 GMT. It was the storms of the 14th in this area which gave the more 
intense short-period falls. Across East Anglia most of the rain fell on the 15th 
starting between 14 and 15 GT and continuing for 8 to 10 hours. This rain 
was generally less intense than that which occurred in Kent and Surrey earlier 
in the day. 


Rainfall on previous days. Rainfall for the first 13 days of September was 
generally less than 50 mm and over quite a large part of the area it was less than 
25mm. Rain fell over the whole area during the first few days, but from the 
5th to gth inclusive most parts were dry. Rain was widespread over the area on 
the roth but amounts were generally low. The 11th was dry over much of the 
area, but a few places in East Anglia recorded between 25 mm and 35 mm. The 
next two days were dry except for parts of Kent and Sussex. As a result of this, 
the soil moisture deficit, which on the 11th ranged from 10 mm to 100 mm, was 
quite large over most of the area. 


Flooding and damage. Serious and widespread flooding occurred over 
Kent, Surrey and parts of Essex, Hertfordshire and Suffolk with localized 
flooding in many other parts of the area. In general, rivers reached their 
peak (maximum height) on the 16th and the water continued to spread on the 
17th causing disruption to all forms of communication and many thousands of 
pounds’ worth of damage. 

The flow in the lower reaches of the River Mole was about four times that of 
January 1955 which was itself the highest recorded between 1930 and 1955. 
In Sussex at Park Mound the River Arun rose 5:2 m in 24 hours, from almost 
zero flow to the highest level ever recorded, and on the Rother the highest level 
previously recorded at Fittleworth was exceeded by 0-3 m. On both of these 
rivers and on many others the significant feature was the extra rainfall on the 
evening of the 15th. By the River Stour in Suffolk the flood was comparable with 
that of 1947. In the Kent River Authority area, flooding in the valleys of the 
Rivers Darent, Cray, Eden, Eden Brook and Bourne was far worse than any 
previously recorded. Besides these valley floods a great deal of flooding occurred 
remote from the rivers and before the rivers overflowed. The flow of a number 
of rivers and streams had to be estimated because the gauging stations were 
either flooded or bypassed by the flood waters. 

The total area of land flooded, excluding the Sussex River Authority area for 
which information is not available, was estimated as 220 km? of which more 
than 100 km? were in the Kent River Authority area and 60 km? in the Thames 
Conservancy area. In this latter region it was the area flooded by just three 
tributaries of the Thames, the Mole, Wey and Loddon. The figure given for the 
Kent River Authority area includes 44 km? of marshland. Again excluding the 
Sussex River Authority area, almost 21 000 properties were flooded. Over 
14.000 of these, of which 10 000 were residential, were in the area of the Thames 
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Conservancy, mainly in the valleys of the Mole and the Wey. It is reported 
that one factory alone suffered a loss of £500 000 including damage and loss of 
production. 

A number of roads were impassable for varying lengths of time including the 
A 23 at Horley and Gatwick, the A 217 at Sidlow, the A 24 at Mickleham, the 
A 3 at the Kingston Bypass, Cobham and the Guildford Bypass, the A 281 at 
Shalford, the A 244 at Hersham, the A 329 near Reading, the A 23 at Brasted, 
the A ai at Darent and the Ag at the confluence of the Rivers Cray and 
Shuttle. At Tonbridge on the 15th all roads northwards were impassable 
owing to local flooding. Several bridges collapsed and others were damaged by 
the volume of water rushing under and round them. More details of damage 
are given in Table II. 

Total cost of damage caused by the exceptionally heavy two-day rainfall of 
14 and 15 September 1968 is impossible even to estimate. The loss of crops on 
agricultural land, and the damage to roads and bridges, and to public and 
personal property cannot even be listed completely. 


TABLE I—HIGHEST RECORDED RAINFALL MEASUREMENTS 


Station National grid 
number reference 


14 September 1968 

294630 Hildenborough TO 51 571488 

295005 Hadlow College TQ 51 628497 

295010 Hadlow Sewage Works TO 51 631491 

297727 Brenchley Gardens, Maidstone TQ 51 759561 

297722 Mill Street, Maidstone TQ 51 761555 

290007 Cross Ness, Erith TQ 51 485809 

237466 Stifford Pumping Station TQ 51 590802 

236974 Marsh Farm Sewage Works, Tilbury TQ 51 654756 

235389 Basildon Sewage Works TQ 51 737908 

235519 Wickford TQ 51 772921 

15 September 1968 

288687 Bromley TQ 51 399693 129°5 

293572 Godstone Sewage Works TQ 51 365505 123°7 

292953 South Godstone Sewage Works TQ 51 363478 124°5 

291253 Priory Gardens, Orpington TQ 51 465666 121°g 

14/15 September (48-h totals) 

237466 Stifford Pumping Station 590802 200°9 

236974 Marsh Farm Sewage Works, Tilbury 654756 201°4 

236632 Canvey Island Pumping Station 800828 155°2 

284974 Earlswood Sewage Works ’ 275483 I51‘l 

293572 Godstone Sewage Works 168-4 

297340 East Farleigh Lock 51 735535 162°3 

297727 Brenchley Gardens, Maidstone TC 759501 157°7 

297347 Barming Water Works, Maidstone q 735548 153°4 

293953 South Godstone Sewage Works k 363478 189°8 
Bromley 399693 190°5 

290200 Sundridge Pumping Station 489556 160°8 

290060 Erith Sewage Works TQ 51 523777 154°4 

290070 Westwood Pumping Station, Surrey TO 51 42554! 150°! 

290241 Chevening Gardens, Sevenoaks TQ 51 484577 162°6 

290412 Stone Street, Kent TO 51 572549 176-2 

294056 Paines Hill, Limpsfield TQ 51 412513 157°8 

294630 Hildenborough TO 51 571488 171-0 

295010 Hadlow Sewage Works TQ 51 631491 165;1 

295945 Yalding, Kent TQ 51 690497 154°4 

297707 Hockers Lane Pumping Station, Maidstone TQ 51 788571 159°6 

297455 College Road, Maidstone TQ 51 757549 165°1 


Note: Many of these values were reported in inches and have been converted to millimetres. 
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TABLE II—FLOOD DAMAGE 


Thames Conservancy area 
Area flooded 


River Mole between Molesey and Horley 
River Wey between Weybridge and Farnham 
River Loddon 


Damage to bridges 


River Mole Betchworth Bridge Loss of both parapets 
Brockham Bridge Loss of one parapet 
Boxhill Bridge Collapse of one arch 
Stoke D’Abernon (rail) Displacement of girders 
Downside Bridge Collapse of two arches 
River Wey Brockford Bridge Collapse of one arch 
Brooklands Partial collapse of span carrying services 
Hoe Stream Godalming (rail) Collapse of one span 
Cranleigh Waters Shalford (A 281) Collapse of arch 


Flooding of properties (provisional figures) 


Local authority* Houses Other premises 
Aldershot M.B. 7 I 

Chertsey U.D. 200 20-30 
Crawley U.D. 150 included 
Dorking and Horley R.D. 200 included 
Easthampstead R.D. 8 o 

Esher U.D. 8000 included 
Farnborough U.D. 24 included 
Farnham U.D. 220 included 
Godalming M.B. 110 10 
Guildford M.B. 300 160 
Guildford R.D. 170 a few 
Hambledon R.D. 400 included 
Hartley Wintney R.D. 24 included 
Leatherhead U.D. 315 50 
Reigate M.B. 200 150 
Walton and Weybridge U.D. 2000-2500 included 
Wokingham U.D. 1000 not known 
Wokingham M.B. 33 2 
Wokingham R.D. 70 included 


Kent River Authority area 
Area flooded 


River Medway 26-9 km? 
River Bourne 2-6 km?® 
River Eden & Eden Brook 15°6 km? 
River Darent 7°2 km? 
River Cray 2:8 km? 
River Shuttle 2°3 km? 
Warden Bay Drainage 0-3 km? 


Marshes 


Allhallows, Grain g°2 km? 
Higham Cliffe 24°8 km? 
East of Gravesend 3°8 km? 
Swanscombe 0-5 km? 
Dartford-Stone 2-3 km? 
Sea-salter 2-8 km? 


Hoo 0-3 km? 





* Local Authority areas are as before 1 April 1974. M.B., U.D. and R.D. denote Municipal 
Borough, Urban District and Rural District. 
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Properties flooded 


Westerham 
Brasted 
Sundridge 
Chipstead 
Dunton Green 
Riverhead 
Sevenoaks 
Kemsing 
Otford 
Shoreham 
Eynsford 
Farmingham 
South Darenth 
Dartford 

St Mary’s Cray 
Crayford 


Edenbridge 
‘ , Hildenborough 
Tonbridge Barden area 


East Peckham 


Yalding 
Maidstone 
Hadlow 


Essex River Authority area 


Are flooded 
Tilbury Marshes 
East and West Tilbury Levels 
Mardyke catchment 


Properties flooded 


Sudbury 

Haverhill 

Keddington 

Sturmer 

Clare and Stoke-by-Clare 
Cavendish 

Long Melford 

Bures 

Cornard 

Stratford St Mary 
Wickford 

Southend 

South Benfleet 

Canvey Island 

Purleigh 
Hawkwell—Hockley—Rochford 
Fobbing 

Rainham 

Loxford Water area 
Mayesbrook area 


Lee Conservancy area 
Area flooded 


Rivers Lee, Mimram, Beane, Rib, Ash, Start 


Properties flooded 


River Beane upstream of Hertford 
River Rib upstream of Hertford 
Bishop Stortford 

River Lee, Hertford to Hoddesdon 


Houses 
110 

18 

12 

61 

10 


some houses 


number not known 

140 

362 

450 

many houses and 
industrial properties 

40 

number not known 


43 


Houses 


Others 
3 


re) 

2 including one bridge 
3 

3 

6 

farm buildings 

farm buildings 

one bridge 

4 + 2 bridges 

2 


2 
° 
258 
re) 


36 


152 


included 


15 


Others 
o 


included 
included 
included 
included 
included 
included 
included 
included 
included 
I 

I 
included 
a few 
industrial premises 
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Greater London Council area 
Area flooded Houses 20°0 km? 
Properties flooded 3000 


Note: Many of the areas given in this table were supplied in acres or square miles and have 
been converted to square kilometres. 


Part II of this paper will be published in October. 
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REVIEWS 


The use of satellite pictures in weather analysis and forecasting, WMO Technical Note 
No. 124, by Ralph K. Anderson, Jerome P. Ashman, Golden R. Farr, 
Edward W. Ferguson, Galina N. Isayeva, Vincent J. Oliver, Frances C. 
Parmenter, Tatiana P. Popova, Rance W. Skidmore, Arthur H. Smith and 
Nikolai F. Veltischev (edited by R. K. Anderson and N. F. Veltischev). 
280 mm X 215 mm, pp. ix + 275, illus., Secretariat of the World Meteoro- 
logical Organization, Geneva, Switzerland, 1973. Price: Sw. Fr. 60. 


This publication is a revision of Technical Note No. 75, produced as a result of a 
recommendation of the Commission for Synoptic Meteorology at its fifth session 
(June-July 1970) that an updated version should be prepared and issued in 
loose-leaf form to facilitate subsequent updating. A welcome feature of this 
new version is the collaboration of U.S.S.R. experts and as a result the incorpora- 
tion of pictures from the cosmos and METEOR satellites and the addition of many 
references to work in the U.S.S.R. 

As might be expected, the opportunity has been taken to present the whole 
subject in a more general and comprehensive way than was possible in 1966 
when Technical Note No. 75 first appeared. Chapter I, ‘General characteristics 
of visible and infra-red images’, serves as an introduction and provides essential 
information on both visual and infra-red images and on the equipment used to 
produce them. Chapter II, ‘Atlas of cloud types, cloud formations and under- 
lying surfaces’, presents the necessary information for the practice of picture 
interpretation. From then on, the emphasis is on the application of satellite 
data, to weather analysis in extratropical latitudes (Chapter III), to synoptic 
analysis in the tropics (Chapter IV) and to the estimation of useful meteoro- 
logical parameters (Chapter V). One of the pleasing new features is the 
arrangement of pictures close to the appropriate text, making for easier 
reference. 

Many forecasters still do not have direct access to satellite pictures but are 
dependent on nephanalyses. Nowhere in this publication is any guidance 
given for the preparation and use of nephanalyses nor indeed any information 
on the preparation and use of coded satellite information. It is well known that 
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satellite information has its maximum value when integrated with more 
conventional data. There is a tendency throughout, particularly in the sections 
on applications to analysis and forecasting in extratropical latitudes to neglect 
this aspect and to treat satellite data as existing in isolation. This is surely a 
mistake, leading inevitably to an over-academic approach. In particular, the 
section on the use of satellite pictures to forecast the evolution of synoptic-scale 
disturbances provides forecasting rules derived solely from picture interpreta- 
tion. Forecasting in extratropical areas is now more dependent on the use of 
numerical models and it is more relevant to know how satellite information 
can be incorporated in analyses to best advantage. Some information on this is 
given in the last chapter. This is obviously a most important subject particularly 
at the present time when the effect of data shortages on the performance of 
numerical models is only too evident. 

Meteorologists everywhere will welcome this revised version and will look 
forward to subsequent extensions incorporating wider concepts of the use of 
satellite pictures in analysis and forecasting. 


T. H. KIRK 


Partial differential equations of mathematical physics, by Tyn Myint-U. 250mm x 
170mm, pp. xiv + 365, illus., Elsevier/Excerpta Medica/North-Holland 
Associated Scientific Publishers, Jan van Galenstraat 335, Amsterdam, The 
Netherlands, 1973. Price: 43.50 Dfl. 


The title is the same as or similar to that of a number of other mathematical 
textbooks that are available and its contents follow fairly closely the traditional 
treatment of the subject. After the introduction and derivation of the equations 
for the physical problems, the canonical forms and general properties of second- 
order partial differential equations are treated and this leads on to the initial 
value or Cauchy problem. After a digression to derive the main results in the 
theory of Fourier series, the method of solution by separation of variables is 
described and illustrated. The Sturm-—Liouville equations and some special 
functions are included in the chapter on eigenvalue problems. Solution of 
boundary-value problems in two dimensions with Dirichlet, Neumann and 
mixed boundary conditions follows. The last three chapters deal with the use 
of Green’s functions to obtain solutions, with integral transforms and with 
equations in three or more space dimensions. 

The contents are selective, and the book cannot rival more comprehensive 
treatises such as Courant Hilbert. But the selection is well made so that the 
argument proceeds without too much need for digression. The approach is 
based on the mathematical problems and their solution, the physical systems 
in which the equations occur entering only briefly in the initial derivation. 
The book is therefore not appropriate for readers who are looking for physical 
insight into the equations or for methods for solving practical problems by 
numerical methods. 

With a text covering such a well-known field the main question that arises 
is probably ‘For what kind of reader can it be recommended ?’. Writing as a 
meteorologist whose mathematics has grown rusty, the reviewer found it 
suitable for the purpose of reminding him of the mathematical basis of many 
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problems arising in meteorology. Almost all the important equations from this 
source are treated, and, from this same point of view, very little of the total 
content could be considered redundant. I suspect most dynamical meteoro- 
logists could read it with considerable profit. For clarity of style and printing 
the text is admirable. Those with the necessary mathematical background will 
find it a pleasure to read. There is a wealth of problems of varying difficulty 
which makes it especially suitable also for students fairly new to the subject. 

While, therefore, it can hardly be recommended as essential reading, this 
book should find a grateful readership among those who are looking for a clear, 
straightforward and comparatively easy account of the mathematical theory 
of the standard partial differential equations which occur in physical problems. 


A. GILCHRIST 


The physics of air-sea interaction, by S. A. Kitaigorodskii. 250mm xX 175 mm, 
pp. vi + 237, illus. (translated from the Russian by the Israel Program for 
Scientific Translations, Jerusalem), John Wiley & Sons Ltd., Baffins Lane, 
Chichester, Sussex, 1974. Price: £7.65. 


Active research in the field of air-sea interaction now involves a range of 
topics too broad to be discussed comprehensively in a single textbook. The 
present monograph is concerned solely with small-scale interactions between 
ocean and atmosphere and its scope is correspondingly limited, but this restric- 
tion in subject matter has given the author enough space to develop some of his 
ideas in considerable detail. Kitaigorodskii has carried out research in a 
number of aspects of small-scale processes near the sea surface and he admits 
that the contents of his book reflect his particular interests; thus, while he has 
written in most cases with considerable authority, there is inevitably selectivity 
also, and some imbalance. 

The text is divided into three main sections. The first occupies roughly half 
the book and is concerned with the dynamics of the near-surface atmosphere 
including vertical turbulence transfers. This part of the book is likely to be of 
most interest to meteorologists. It begins with an introductory chapter on 
flow over a solid boundary and follows with descriptions of the changes which 
result in the flow when the rigid boundary is replaced by the waving sea surfaces. 
The section concludes with empirical data on the exchange coefficients for 
vertical transfers of momentum, heat and water vapour, and the effects of 
departures from neutral density stratification. Of particular interest is the dis- 
cussion of the complex (and as yet undiscovered) relation between drag 
coefficient, wind speed and sea state. 

The second main section is concerned with ocean waves. Dimensional 
arguments are used to derive the principal features of spectra of fully developed 
waves, and of waves after limited fetch or limited duration of wind. A brief 
description is given of the Miles—Phillips theory of wave development. 

The final section deals with vertical mixing processes in the upper ocean. 
This is probably the least satisfactory part of the book, and Kitaigorodskii 
himself writes “The . . . section appears to be premature’. There is a consider- 
able amount of empiricism and while several problems are formulated, few 
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satisfactory answers emerge, even, for example, when important factors such as 
existing stratification or time dependence are omitted from discussions of wind 
mixing. 

The translation suffers from a few flaws which careful editing by a micro- 
meteorologist or oceanographer would have eliminated. For example, it is 
unusual to talk of a ‘deep minimum’ in the power spectrum of wind velocity 
(p. 24). The retention of conventional Russian notation hinders easy reading: 
examples are P for air density and gg, Wx for the vertical fluxes of heat and 
water vapour. However, in spite of the volume’s restricted scope and its 
blemishes it is a welcome addition to the literature of air-sea interaction. The 
price is no higher than one expects in these inflationary times. 


N. THOMPSON 


OFFICIAL PUBLICATION 


The following publication has recently been issued: Annual report on the Meteoro- 
logical Office 1973, London, HMSO, 1974. Price: £1.10. 


The Director-General, in the foreword to his report for the year ending 
31 December 1973, highlights the marked increase in demand for meteorological 
information from industry and the general public. The number of direct 
inquiries to the Office reached an all-time record of 1.75 million, while another 
all-time record of 15.5 million calls were made to the Post Office Automatic 
Telephone Weather Service. 

The computer-based forecasts were further refined during the year, but the 
Director-General stresses the difficulties of maintaining an adequate network 
of observations, particularly over the oceans, and warns that any deterioration 
of the network could have an adverse effect on the accuracy of weather forecasts. 

Activity in research was maintained at a high level. Work on the develop- 
ment of mathematical models to simulate and perhaps predict the behaviour 
of the atmosphere on relatively long time-scales continued. Anomalies in sea 
surface temperatures can have far-reaching effects on the development of 
weather patterns and one very interesting experiment strongly suggested that 
an area of abnormally warm ocean off West Africa was an important factor in 
producing the very cold winter of 1962-63. 

The production of nitrogen oxides in the stratosphere by supersonic transport 
aircraft and their possible effects on ozone concentration at high levels have 
aroused considerable controversy. Calculations by Meteorological Office 
scientists indicate that the quantities of nitrogen oxides injected into the strato- 
sphere by nuclear test explosions in the 1950s and 1960s were equivalent to those 
likely to be produced by 1000 Concordes each flying 7 hours a day. Moreover, 
analysis of world-wide ozone reports before and after the nuclear explosions 
fails to reveal any detectable changes in ozone content. 
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HONOUR 


The following appointment to the Imperial Service Order was announced in 
the Queen’s Birthday Honours List, 1974: 
Mr H. C. Shellard, Principal Scientific Officer, Met O 3. 


OBITUARY 


It is with regret that we have to record the death of Mr P. Jenkins, Assistant 
Scientific Officer, Met O 2, on 31 May, 1974. 
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